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Abstract 



C^ ■ In this paper we study elliptic curves which have a number of points 

whose coordinates are in arithmetic progression. We first motivate 
this diophantine problem, prove some results, provide a number of 
interesting examples and, finally point out open questions which focus 
on the most interesting aspects of the problem for us. 
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1 Introduction 



During this paper, all standard results unless otherwise explicitely 
stated are taken from fH". Extensive use has been made of Maple V + 

^>; ■ APECS (by Ian Connell, McGill University). 

^ . We will deal with elliptic curves defined over a field K hy a Weier- 

H ' strass equation, that is 

E:Y^ + aiXY + a-sY = X'^ + osA^ + a^X + as, a^ G K. 

We win denote, as usual, E{K) the locus of the above equation, 
together with the point at infinity, O = (0 : 1 : 0). 

Changes of variables preserving this form are those given by 

X' = u^X + r, Y' = u^Y + sX + t; 



and we will consider that two equations related by such a change of 
variables represent the same curve (equivalently, we will deal with 
elliptic curves up to so-called Weierstrass changes of variables) . 

Consider then Po)---)-fn G E{K), with Pi = {xi,yi) such that 
xq, ..., x„ is an arithmetic progression. Then we say that Pq) •••> -fn are 
in x-arithmetic progression (x-a.p.) and also E is said to have an 
x-arithnietic progression of length n+ 1. From the previous remarks, 
this does not depend on the Weierstrass equation considered. 

The same definition goes for ^-arithmetic progressions (y-a.p.). 
However, in this case, changes of variables (even those which preserve 
Weierstrass equations) can create and destroy y-arithmetic progres- 
sions. 

Example.— Let us consider the equation over Q 

E-.Y'^ - —XY + —Y = X^ X^, 

16 64 64 ' 



which verifies that 



1 -4 



-1 -3\ / 5 -2\ /I -1 



8' 1287' V 32 '128/' V64' 128/ ' V 32' 128 

(ei'^)' (ei'm)' (le' m) "^^^^^^ 

The reader can easily check that, after the change of variables 

Y' = Y + X, 

for instance, the corresponding points are not in y-a.p. Hence we can 
properly talk of x-a.p. in a curve, but if we speak of y-a.p. in a curve 
we must bear in mind that we are considering a specific equation. 

This paper studies with elliptic curves which have a simultaneous 
arithmetic progression. First we need a proper definition of these 
progressions. Let us consider PQ,...,Pn as above. If we ask for both 
xo, ..., x„ and yo, ...,yn tohe arithmetic progressions then the problem 
is far too easy, as Pq, ...,Pn must be collinear and hence n < 2. Lots 
of examples can be found with this property; for instance all curves 
in the family 

E{b) : y2 + (25 _ l)XY + bY = X^ - bX"^ 
have the arithmetic progression (0, —6), (6,0), {2b, b). 



Definition.— With the above notation, Pq, ...,Pn are a simultaneous 
arithmetic progression (s.a.p.) or the equation E is said to have a 
simultaneous arithmetic progression if: 

(a) xq, ...,Xn are in arithmetic progression (called the support of the 
s.a.p.). 

(b) There exists a permutation a in the symmetric group of n + 1 
elements Sn+i such that ya{o), ■■■lUa^n) ^^^ in arithmetic progres- 
sion. 

The definition is clearly symmetric: it is equivalent (up to point or- 
dering) to saying that yo, ..., yn are in arithmetic progression and there 
exists fi € Sn+i such that x^(-o), ...,x^(-„) are in arithmetic progression, 
but this version showed more useful for computational purposes. 

With this definition, at least three straight problems arise: 

(a) The detection problem: Given an elliptic curve, does there exist 
an algorithm for deciding whether it contains or not a s.a.p. of 
length n (giving as well a change of variables if needed)? 

(b) The subsequence problem: If an elliptic curve has a s.a.p. of 
length n + 1, does it possess a s.a.p. of length n? (Note this is 
not at all trivial from the definition). 

(c) The bound problem: Is there a bound for the possible lengths of 
s.a.p. in elliptic curves? 

Trying to understand these three problems, we have developed 
some computational methods (actually, two partial answer to the de- 
tection problem) whose application may shed some light. Nevertheless 
the results achieved can be considered only as a first step towards a 
fully satisfactory understanding of these sequences. We have managed 
to prove the following results: 

Theorem 1.— Given an elliptic curve with an x-a.p., there exists an 
algorithm which decides whether or not the curve also has a s.a.p. 
with the given x-a.p. as support. 

Theorem 2.— There are integers n such that there are examples of 
s.a.p. of length n which do not contain any s.a.p. of length n — 1. 

Theorem 3.— There are no elliptic curves defined over Q with s.a.p. 
of length 7. There are only finitely many non-isomorphic curves de- 
fined over Q with s.a.p. of length 6. 



We will finish this introduction with a word on motivation. At 
first our interest was drawn to this subject by the articles of Bremner- 
Silverman-Tzanakis (|4|) and Bremner (|3 ). Apparently these papers 
had their starting point in the featuring of x-a.p. as by-product of a 
latin square problem (see more on this in j^H]). However, highly in- 
teresting results were sketched in both papers around the relationship 
between the existence of arithmetic progressions on a certain elliptic 
curve and its rank. In this same line a paper by Campbell appeared 
(0) pointing out far-reaching questions, probably too difficult for the 
state-of-the-art . 

The history of the problem, though, can be traced back quite fur- 
ther, as (for the specific case of Mordell curves) it was treated previ- 
ously by S.P. Mohanty (8) who studied x and y-a.p. separatedly and 
by Lee and Velez ([2|) who first treated s.a.p., if only in the naive form 
mentioned above, without permutations involved. The motivation for 
these first works was, as many other times in the history of number 
theory, purely diophantine. 

We became interested in this specific problem while trying to im- 
prove Bremner's record of longest x-a.p. by narrowing the search^. 
Our first attempts were shown in [S], using a specific kind of s.a.p., 
which allowed us to find examples of s.a.p. of length 5. These methods 
were not at all exhaustive, as it accurately pointed out by Bremner 
in his MathSciNet review. After this work, we feel that some of the 
problems posed are worth a closer look and the setup remains chal- 
lenging. As Bremner points out in 3 : "Questions in number theory 
that interrelate two group structures are easily posed, but often lead to 
intractable problems". 

2 The detection problem 

Let us consider a set of points Pq = {xo,yo), ..., Pn = (xmlln) in an 
elliptic curve, defined over i^ by a Weierstrass equation: 

E ■.Y'^ + aiXY + a^Y = X^ + 02^^ + 04^ + a^, with Oj G K. 

Let us suppose the points Pq, ..., Pn G E{K) to be in x-a.p. We are 
interested then on finding, if there exists any, a change of variables, 
preserving the Weierstrass form of E, which transforms Pi = {xi,yi) 



^The longest x-a.p. found in an elliptic curve has 8 terms (^); for y-a.p. the record 

is 7 so far m). 



into Pi = {x[,y[) so that Pq,...,PI^ is a s.a.p. on the corresponding 
equation E' . This change of variables must be of the form 

X' = u^X + s, Y' = u^Y + rX + t. 

If Xj = a + i ■ d, and we want y'^(Q\, ■•■^y'a(n) ^° ^^ ^'^ arithmetic 
progression for some a € 5*71+1, then it must hold 

x[ = u'^{a + id) + s 

y[ = h + a{i)d' , for i = 0, ...,n. 

We can take, with no loss of generality, it = l, s = 0, t = 0. This 
involves only choosing an appropriate reference system by translation 
and scaling (which would not affect s.a.p. in any case). Then we have 



X^ — Xi 



a + id 



y'i = Vi — r{a + id) =h + a{i)d' , for i = 0, ..., n. 

for some b, d' € K. 

These last identities can be written as a system of n + 1 linear 
equations in r, 6 and d' , with matrix 



A* 



f a + Od 1 a{0) 
a + ld 1 cj(l) 

y a + nd 1 cr{n) 



Vi 
Vn j 



Note that the y-sequence yo,...,yn is not an arithmetic progres- 
sion if and only if first, second and fourth columns are independent; 
equivalently 



3s E {2, ...,n} such that 






1 


yo 


1 


1 


yi 


s 


1 


ys 



y^O. 



Algorithm 1.— Our first detection algorithm is based on the fact 
that, the existence of a solution to our system (that is, the existence 
of a s.a.p.) is equivalent to A* having rank 3. The formal algorithm 
goes like this: 

Input Data: E, xq, ...,Xn (equivalently E, xq, n, d). 

Step 0: Choose a suitable set {yo,...,yn} such that {xi,y,i) G E{K). 



Step 1: (Fool-proof checking) Check whether {t/Q; ■■■iVn} is an arith- 
metic progression. If so, we are finished; if not, find an s, 2 < s < n 
as above. 

Step 2: For any a £ Sn and any i £ {2,...,n}, i j^ s, compute the 
minor formed by the first, second, s-th and i-th rows of A*. 

Step 3: If, for some a £ Sn the n — 2 minors are null, solve the system 
to find r, b and d'. If not, back to step 0. 

The main inconvenience of this procedure is its needing of 2"^^(n+ 
l)!(n — 2) determinant computations, as there are (n + 1)! possibilities 
for a and 2 possibilities for each yj. So we will try to find a more 
efficient procedure, although this set-up will prove useful later on. 

Algorithm 2.— Consider the affine points Qq = (0, yo, o"(0)), ..., Qn = 
{nd,yn,cr{n)) G A^(K). Note that, what we need from all them, in 
order to have a s.a.p., is to be in the same plane. This is the basis 
for our algorithm, whose input data and steps and 1 are identical to 
the previous one: 

Step 2: For each {i,j, k} C {0, 1, ...,n}, we consider the plane 

■^ijk = {{0,yo,'i),{d,yi,j),{sd,ys,k)), 

Step 3: For / = 2, ...,n and I ^ s, we intersect the line 

x = ld, y = yi 

with the plane vtjjA;- 

Step 4a: If any of these intersections gives a point {ld,yi,zi) such 
that zi ^ {0, 1, ..., n} or z/ is equal to another zi', then {zq, ■■■,Zn} does 
not correspond to {o"(0), ...,a{n)} for any a G Sn+i- Back to step 2, 
change the plane and repeat the process or back to step if all planes 
have been exhausted. 

Step 4b: If we find out a set of points Qo = {xo,yo,zo), ...,Qn = 
{xn-, Un, Zn) with Zj = a{i) for i = 0, ..., n and a € S'n+i, then a allows 
us to have a solution r, b, d' of our system. 

As for computations is concerned, note that we have (n + l)n{n — 
l)/6 possibilities for TTijk, and for each plane we have, at most, n — 2 
intersections. This, together with the 2"+^ possibilities for each yj, 
means a saving of around (n — 3)! computations. 



The implementation of both algorithms shows the time difference 
is not huge (as n dos not go very far), but it already grows significantly 
for n < 5. 

From now on, we will note a € Sn+i by a = (ao...a„), meaning 
o-(O) = ao,...,cr(n) = a„. 

Example.- Let the curve Y'^ = X"^ - 112X + 400, defined over Q, 
which has the following x-arithmetic progression of length 4: 

Xk : -4,0,4,8, 

and there are 4 of the 2^ y-sequences that lead to simultaneous arith- 
metic progressions. 



yo 


yi 


y2 


ys 


a 


y'o 


y'l 


y'2 


y'3 


28 


-20 


4 


4 


(1023) 


4 


-20 


28 


52 


-28 


-20 


-4 


4 


(0213) 
(1032) 
(1302) 


-44/3 
-18 

-84/5 


-20 
-20 
-20 


-52/3 

-14 

-76/5 


-68/3 

-16 

-92/5 


28 


20 


4 


-4 


(0213) 
(1032) 
(1302) 


44/3 

18 

84/5 


20 
20 
20 


52/3 

14 

76/5 


68/3 

16 

92/5 


-28 


20 


-4 


-4 


(1023) 


-4 


20 


-28 


-52 



Now this x-arithmetic progression can be extended to one of length 



5: 



Xfc: -4,0,4,8,12, 

and there are 2 of the 2^ possible y-sequences that lead to simultane- 
ous arithmetic progressions. 



yo 


yi 


y2 


ys 


y4 


a 


y'o 


y'l 


y'2 


y's 


y'4 


-28 


-20 


-4 


4 


28 


(13240) 


-44/3 


-20 


-52/3 


-68/3 


-12 


28 


20 


4 


-4 


-28 


(13240) 


44/3 


20 


52/3 


68/3 


12 



The equation for both cases is 



n 20 o 

y2 XY = X^ 



—X'^ - \V1X + 400. 



Now, if we try to repeat the procedure for length 6 with 
Xfc: -4,0,4,8,12,16, 



we find that none of the 2^ possible y-sequences leads to a simultane- 
ous arithmetic progression. 

Open problem 1: Find a procedure for deciding whether an elliptic 
curve has an x-a.p. of given length. 

The most interesting results on this line are the parametrizations 
by Bremner in [^j which will be used afterwards in this paper. However 
they are still far from being useful from a computational point of view. 

3 The subsequence problem 

The programs developed in the previous section (specially the second 
one) were of great help with testing the examples we had created 
with the techniques shown on 6^ and also with creating new ones. 
The counterexamples announced in Theorem 2 were product of these 
extensive calculations. Here we present the simplest one. 

(Counter)Example.— Consider the following elliptic curve over Q, 
in Tate normal form, 

/25 -2\ o 25 2 q 2 o 

El—, — : y2 + _XY --Y = X^ + -X^, 
V21' 7 y 21 7 7 ' 

which has the x-arithmetic progression of length 5: 

• Z^ Zi Z^ n ^ 

Using the above procedure we find a y-sequence that gives simul- 
taneous arithmetic progression: 

4 16 92 4 
^'^^ 7' 147' 147'°' 21 • 

There is only one permutation a which passes Algorithm 2 and, 
henceforth, allows the change of variables, a = (20413). The y'- 
arithmetic progression is 



, 8 -8 24 16 




^fc- 49' 49 '49' '49' 




for the equation 




^^+21^^-7^ = ^'+147^- 


" 147^- 



This way, we have found a simultaneous arithmetic progression 
of length 5 not containing a simultaneous arithmetic progression of 
length 4, because the permutation involved is not the extension of an 
54 permutation. In our many calculations these are singular cases: 
first of all, the permutation found is seldom unique and, among the 
collected ones, we usually find an extension of some ^4 permutation. 
But the fact is that these counterexamples happen, which, by the way, 
carry an additional difficulty for all arguments involving induction. 

Interestingly enough, there are other suitable s.a.p. of length 5 
with the same support (that is, other choices for the yi) which happen 
to have subsequences of length 4. 

Open problem 2: Given an elliptic curve with a s.a.p. of length n, 
is there always a s.a.p. of length n — 1 whose support is contained in 
the support of the given one? 

We have found no examples to support a negative answer to this 
question which, by the way, may serve as a weak induction result. 

4 The bound problem 

Our final look will be to the bound problem. As it was pointed out by 
Bremner in [3j, this kind of problems tend to become unmanageable 
quite quickly. From our many attempts, we will describe here the most 
successful of them all; which relies on a parametrization of curves with 
x-a.p. due to Bremner (3;) (here slightly changed for our purposes). 
In what follows we will assume K = Q. Note that all the previous 
arguments do not rely on the base field at all. 

First of all, we will parametrize elliptic curves in short Weierstrass 
form 

Y^ = X^ + AX + B 

with four points in x-a.p.; 

Po = {a,yo), Pi = {a + d,yi), P2 = {a + 2d,y2), P3 = {a + 3d,y3). 

Now, we consider the four polynomials Fq, ...jF^ given by 

Fi = yf - (a + idf - A{a + id) - B, 

in Cl[yQ,yi,y2,y3,a,d, A, B], and compute a Grobner basis of the ideal 
(Fo, ..., -F3). The tdeg ordering in Maple V gives a basis of ten ele- 
ments which can be used for computing a,A,B taking d,yQ, ...,y3 as 



9 



parameters: 
^ = ^ (^0 - ^^0^1 + 6^22/0 + ylvl + 21yi - 39yi2/? + 6yiy?+ 

.4 n„,2„,2 , „,4\ D/QCj4 



B = ■^^{ytvl + ^yiyl + yiyl-'^yly'-^ylylyl + '^^ylyf 

-^ylvhl - Uyjylyl + yfy^ + y^^ - dy^yf + 20yf 
-21yfyj + 4y^y2 + 20^2' - 21y?y| + 24y22/2' 
-I2ylylyl) = Q/Q'<f 

a = ^ (-2%' + %? - %! + yt) = -^/6rf' 

and, in addition, the first member of the basis is 

-yl + Qd' + yl-2,yl + ^yl 

Now, making the scaling with u = Qd we obtain the following 
parametrization: 

Curve: Y"^ = X^ - Q^P + 6^Q 
First term: 3Gd?a = —6R 
Difference: 36d^ = 6y| - 6y§ + I8yf - 18y^ 

We will use from now on A,B,a,d for these new polynomials. It 
is interesting noting that the points in x-a.p. are now 

Pi = {a + id, ±Qyid), for i = 0, ..., 3. 

We will try to produce curves with a s.a.p. of given length with a 
variant of Algorithm 1 which we will illustrate with the case of length 
6. In fact, using this procedure we might compute all curves with such 
a s.a.p. in contrast with the lack of exhaustiveness of [S]. If we want 
points P4 = (a + 4d, Z4) and P5 = (a + 5d, Z5) to be in the curve it 
must hold 



n) 



Z4 = ±m^4.yl + Ayj-y^o-6y^{-yl + y^o-3yi + 3y. 

Z5 = ±36v'-4y2 _ 20yl + Ibyf + lOyU-yl + ^o " 3y? + 3yi) 

Hence, for the sake of consistency, we will call 

yj = ^yl + ^yl-yl-Qyl 

yl = -4y2 - 20yi + 15y? + lOyi 
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M 



and our new points will then be Pj = (a + id, ±6yid) for i = 4, 5. 

Example.— We will show how to proceed using the 3-cycle a = 
(210345). After the above remarks, we may use the matrix M, given 

by 

/o 1 

1 1 

2 1 

3 1 

4 1 
V5 1 

instead of the original matrix A* from Algorithm 1; and ask M to 
have rank 3, as we are assuming d j^ 0. Note that {i/q, ...,^5} are not 
the y-coordinates of the points Pq, ...jP^. As all the minors are linear 
polynomials on yo,...,y5 and we have also two quadratic relations, it 
is not surprising that the complete solutions are two linear varieties, 
actually a plane and a line, given by the following parametrizations: 



1 


yo \ 


2 


yi 





y2 


3 


ys 


4 


yA 


5 


yb J 



yo, 



2/4 + 3yo y4 + yo 3y4 + yo 



y4, 



72/0 ._ 



m 

19' 



-15yo -27yo 



5y4 - yo 
4 

-39yo\ 

19 r 



19 19 19 19 

from which the first one only contains points inducing d = and 
therefore must be discarded. In fact, these trivial solutions appear in 
all cases, which is clearly a by-product of our previous assumptions. 
Now we make the substitutions induced by the second parametriza- 
tion, obtaining 



A 



-7840512 



-yo' 



B 



8449090560 



-yo 



1536 



130321 "" 47045881 ^ 361 

and the linear system given by A* has solution 

_ -60yo , _ 98208y3 



d' 



■yo, 



576y3 



d 



48 
361 



yo> 



19 ' 193 ' 193 ■ 

This gives, after the corresponding substitution, the equation 



Y' 



UOyo 



XY = X" 



3600y^_^2 



19 361 

which has the following s.a.p. 



7840512y^ 8449090560yg 

-X + 



130321 



47045881 



/l536yg 97632y|]\ /l584y§ 97056y|]\ /l632yg 978208yg^ 



192 



193 



192 



193 



192 



193 
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(l680yl 96480y|\ /l728y| 95904y|\ /l776y| 95328y|\ 
1^ 192 ' 193 J' [ 192 ' 193 J' [ 192 ' 193 j ' 

All these curves are isomorphic to (an easy) one given by the case 
yo = 19/2; 

y2 _ QOYX = X^ - 900^2 - 490032X + 132017040, 

having the sequence {(384, 12204), (396,12132), (408,12276), (420,12060), 
(432,11988), (444,11916)}. 

We have not computed all curves with s.a.p. of length 6, although 
we have bounded the number of curves by 19200 cases, using the 
previous computations with all possible sign and permutation choices, 
counting only the number of possible solutions, that is, cases where 
the line does not induce d = 0. 

To be precise, only half of the sign choices have to be considered, 
as every arithmetic progression of difference d is also an arithemtic 
progression of difference —d, and hence, every curve appears at least 
twice, for a pair of inverse choice of signs and permutations. 

Even so, not all of these cases are non-isomorphic elliptic curves; 
there might be isomorphic curves among them as well as genus 
curves. In the appendix we have given some explicit data for the first 
100 curves actually computed with this method where repeated curves 
already appear (in fact, there are only 56 non-isomorphic curves). 

As a side remark, the distribution of the possible is extremely reg- 
ular: there are 600 allowed permutations (that is 600 lines not giving 
d = 0) for every sign choice and only sign changes were allowed most 
for the solutions for a fixed permutation. By the way, these differences 
ususally disappeared when finding the solutions to the system given 
by A*. This, together with the repeated cases shown in the appendix, 
gives a heuristic estimation of only around 350 non-isomorphic curves, 
but filling the details of such a list is beyond our computational pos- 
sibilities so far. 

As for length 7 is concerned our procedure shows there are no 
solutions, as all induce d = 0. This case exhausts the possibilities 
of computer checking, at least with these methods, as it took around 
20 hours of CPU (which implied four days in real time) and, more 
constraining, 2^^ Mb of stack memory. For considering this attack to 
length 8, these figures should be multiplied at least by 16 (8 for the 
number of permutations and there are twice as many minors now) , let 
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alone the additional difficulty of adding a new quadratic polynomial 
to the system, which is not easy to measure. 

Open problem 3: Find a universal bound for the length of a s.a.p. 
on elliptic curves over Q. 

Note that an affirmative answer to the open problem 2 would mean 
6 is the answer to open problem 3. 

Appendix: Examples of curves with s.a.p. 
of length 6 

Please note that this stream was computer-generated taking as a 
unique choice the permutation. Due to this, some of the curves (for 
instance examples 009 and 010) appear more than once, as they have 
different s.a.p. Also one may find isomorphic curves (as 001 and 002). 
As we said above, only 56 non-isomorphic curves can be found in this 
table, but we have preferred to leave as it came, as we feel it illustrated 
better the phenomenon. The entries of the table after each equation 
are: 

1) Permutation (a): noted as above by (cj(0)cj(1)...o"(5)). 

2) Numerical data (N.D.): The set [a,d,b,d'] which fits the equa- 
tion. 

3) Rank (r): The rank of the curve, computed with APECS (< 
means APECS failed to actually compute the rank, in which case the 
best upper bound given is shown). All curves have trivial torsion 
group. 

001 Y^ - 180yX + SIOOX^ -X^ + 4892251392X - 134063884477440 
a = (321450), N.D. = [66432, -20304, 13044672, -5725728], r = 5 

002 Y^ + 180YX + 8100X2 - X^ + 4892251392X - 134063884477440 
a = (054123), N.D. = [-35088,20304, -13044672, 5725728], r = 5 

003 Y^ + 20YX + lOOX^ - X^ + 36478512X - 82321246080 

a = (423150), N.D. = [5724, -1584, -367704, 139392], r = 5 

004 y2 - mYX + 900X2 - X^ + 86832X - 8864640 
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a = (534201), N.D. = [-324, 144, 9288, -3456], r = 2 



005 Y^ - 20yX + lOOX^ - X^ + 466992X - 549797760 



a = (045312), N.D. = [1308, -432, -33576, 3456], r = 3 



006 y^ - 2YX + X^ - X^ + 238707X - 41709006 

a = (150423), N.D. = [-513, 180, 8487, -3600], r = 3 



007 y^ - ISOyX + 8100^^ - X^ + 51432192X - 368371860480 
a = (034521), N.D. = [-9312, 4752, -638496, 313632], r = 4 

008 y^ - 612yX + 93636X2 - X^ + 26962612992X - 1882111797863424 
a = (145032), N.D. = [-164832, 66960, -19904832, 12454560], r = 5 

009 y^ + 396yX + 39204X2 - X^ + 388106595072X - 88686989876929536 
<T = (250143), N.D. = [-715488,235440,459510624, -153977760], r = 5 

010 y2 - 576yX + 82944X2 - X^ + 271059091200X - 49072046238950400 
o- = (134502), N.D. = [734160, -265680, -415035360, 130714560], r = 5 



a = (350124), N.D. = [-594240, 265680, 238537440, -130714560], r = 5 

y2 + 396yX + 39204X2 - X^ + 388106595072X - 88686989876929536 
a = (305412), N.D. = [-715488, 235440, -310378176, 153977760], r = 5 

y2 - 612yX + 93636X2 _ x^ + 26962612992X - 1882111797863424 
a = (410523), N.D. = [-164832,66960,42367968, -12454560], r = 5 



014 y2 - 180yX + 8100X2 - X^ + 51432192X - 368371860480 
a = (521034), N.D. = [-9312,4752,929664, -313632], r = 4 



015 y2 - 3ooyx + 22500x2 - x^ + 894074112X - 19561912750080 

a = (230145), N.D. = [-30384, 11808, -6045408, 2904768], r = 6 

016 y2 - 60yX + 900X2 _x^ + 48O45312X - 189087436800 

a = (341250), N.D. = [-7344, 3744, -714528, 292032], r = 5 



017 y2 - 60yX + 900X2 _ x'-^ + 48O45312X - 189087436800 

a = (503412), N.D. = [11376, -3744, 745632, -292032], r = 5 
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018 Y^ + 300YX + 22500X^ - X^ + 894074112X - 19561912750080 
o- = (014523), N.D. = [28656, -11808, -8478432, 2904768], r = 6 

019 Y^ - 1800YX + 810000X2 - X^ + 7410269014272X - 8108956611489899520 
a = (302145), N.D. = [2763696, -577584,3502719072, -1323822528], r < 6 

020 Y^ + 2880YX + 2073600^^ - X'-^ + 403246536682752X - 3207719147336296058880 
o- = (413250), N.D. = [-829968,5312736,62302243392, -30155089536], r < 6 

021 Y^ + 4680YX + 5475600^^ - X^ + 3994541989632X - 6207305561351930880 
a = (524301), N.D. = [-2327712, 728784, -3390388704, 2107643328], r < 7 

022 Y^ - 2880YX + 2073600X2 - X'-^ + 8792422324992X - 14409886139859502080 
a = (035412), N.D. = [7563792, -2228688, -8618074272, 1791865152], r < 11 



a = (140523), N.D. = [-42040752, 12290400, 265952884032, -83918851200], r < 7 



024 Y^ - 6912YX + 11943936X2 - X^ + 214026269587200X - 1973672371225801958400 

a = (250134), N.D. = [-9752640, 7294320, 76016037600, -26872274880], r < 8 



a = (532140), N.D. = [8064, -3168, 1169856, -418176], r = 4 



a = (043251), N.D. = [11544048, -3013920, -34095021120, 10813944960], r = 4 

027 y2 + 2616YX + 1710864X2 - X^ + 3327107104512X - 3360005144504534016 
o- = (154302), N.D. = [-515808,403920, -2454087456,906396480], r < 6 

028 y2 - 3072YX + 2359296X2 - X^ + 2328753038487552X - 49786655684722942869504 
a = (205413), N.D. = [-55130112, 19313280, 240362650368, -111012733440], r < 6 

029 y2 - 4944YX + 6110784X2 - X^ + 14552162687232X - 21391150302252205056 
(7 = (310524), N.D. = [1164672,340560, 5468006304, -44953920], r<7 

030 y2 - 960YX + 230400X2 - X^ + 19321742592X - 1033886822584320 
o- = (421035), N.D. = [60912, 6048, 38636352, -72576], r = 5 
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031 Y^ + mYX + 400X2 - X^ + 8854272X - 8259978240 

a = (153420), N.D. = [-2448, 1584, 242208, -69696], r = 4 



a = (204531), N.D. = [-1152, 720, -44064, 14400], r = 3 



033 Y^ - 64YX + 1024^^ - X^ + 1389312X - 394896384 

a = (420153), N.D. = [2448, -720,27936, -14400], r = 3 



034 y-^ - 40yX + 400X-^ - X"' + 8854272X - 8259978240 

a = (531204), N.D. = [5472, -1584, 106272, -69696], r = 4 

035 Y^ - lOSOyX + 291600X2 - X^ + 4379037151491072X - 136063103467431710822400 
a = (240315), N.D. = [-77988192,33054048, -438167226816,226486336896], r < 8 

036 y2 _ 720YX + 129600X2 - X^ + 9969629720832X - 13778174775900128256 
a = (351420), N.D. = [-3640128, 1648080, -6422790240, 2155688640], r = 1 

037 y2 + 5544yX + 7683984X2 _ x^ + 422337203926272X - 3401406251614088487936 
a = (402531), N.D. = [-18956688, 10568880, 102961496736, -29550588480], r < 7 

038 y2 _,_ losorx + 291600x2 - X^ + 4379037151491072X - 136063103467431710822400 
a = (513042), N.D. = [87282048, -33054048,438167226816, -226486336896], r < 8 

039 y2 - 720YX + 129600X2 - X^ + 9969629720832X - 13778174775900128256 
a = (024153), N.D. = [4600272, -1648080, -6422790240, 2155688640], r = 1 

040 y2 - 55UYX + 7683984X2 _ x'^ + 422337203926272X - 3401406251614088487936 
a = (135204), N.D. = [33887712, -10568880, -102961496736,29550588480], r < 7 

041 y2 + 20yX + 100X2 - X^ + 466992X - 549797760 



o- = (342015), N.D. = [-852, 432, 16296, 3456], r = 3 



042 y2 - 60YX + 900X2 _ x^ + 86832X - 8864640 

a = (453120), N.D. = [396, -144, -7992,3456], r = 2 

043 y2 - 20yX + 100X2 - X^ + 36478512X - 82321246080 

a = (504231), N.D. = [-2196, 1584, -329256, 139392], r = 5 

044 y2 - 4yX + 4X2 _x-^ + 3819312X - 2669376384 
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a = (231504), N.D. = [1548, -720, -76104, 28800], r = 3 



a = (402315), N.D. = [14071152, -3627936,35652163392, -14061879936], r < 7 

046 y^ + 360yX + 32400X2 - X^ + 27909792000X - 1526822144640000 
a = (513420), N.D. = [-99600, 75600, 73720800, -22680000], r < 5 

047 Y^ - 7056yX + 12446784^^ - X^ + 129437671756032X - 543091962160151110656 
a = (135042), N.D. = [17527152, -4944240, -5119105824, 6467065920], r < 5 

048 Y^ - 7488yX + 14017536X2 - X^ + 841545692773632X - 15278101189303210758144 
a = (240153), N.D. = [-34464912, 10864800, -202374030528, 65579932800], r < 8 



049 Y^ - 6696yX + 11209104X2 - X"' + 1767944739520512X - 27067552170879621206016 
a = (351204), N.D. = [-46912128, 22001760, 291906110016, -122241778560], r < 7 

050 y2 - 90yX + 2025X2 - X^ + 3434834187X - 77407635019590 



a = (053421), N.D. = [39747, -2376, 3706047, -627264], r = 5 



a = (104532), N.D. = [-21237, 11880, 3783861, -1425600], r = 4 



a = (320154), N.D. = [152652, -47520, 26753112, -11404800], r = 4 



a = (431205), N.D. = [27867, 2376, 569727, 627264], r = 5 



054 y2 + 396YX + 39204X2 - X^ + 388106595072X - 88686989876929536 

o- = (214503), N.D. = [461712, -235440, -310378176, 153977760], r = 5 



055 y2 - 612yX + 93636X2 _ x^ ^ 26962612992X - 1882111797863424 
a = (325014), N.D. = [169968, -66960,42367968, -12454560], r = 5 

056 y2 _ 180YX + 8100X2 - X^ + 51432192X - 368371860480 
a = (430125), N.D. = [14448, -4752, 929664, -313632], r = 4 

057 y2 - 48yX + 576X2 _ X'^ + 9262512X - 8383430016 

a = (250341), N.D. = [-3348, 1440, -123552,43200], r = 2 

17 



a = (412503), N.D. = [963, -360, 11556, -5400], r = 2 



059 Y^ - QOYX + 900X2 - X^ + 767232X - 253808640 

a = (523014), N.D. = [1008, -288, 16416, -1728], r = 3 

060 y2 - 6oyx + 900^^ - x^ + 767232X - 2538O8640 

a = (145230), N.D. = [-432, 288, 7776, 1728], r = 3 

061 y2 - 5184yX + 6718464X2 - X^ + 142783615579392X - 708085096419827994624 
o- = (352041), N.D. = [13189008, -5387040, -40963299264, 13898563200], r = 4 

062 y2 - 2142yX + 1147041X2 - X^ + 15808640228307X - 39293749322977053294 
o- = (403152), N.D. = [5885823, -1961820, 10341141993, -4072738320], r < 7 



a = (514203), N.D. = [-243780, 110160, 98574840, -44945280], r < 5 

064 y2 - 396YX + 39204X2 - X^ + 16655166000X - 690091901769600 

a = (025314), N.D. = [250140, -71280, -60813720, 18817920], r = A 



a = (130425), N.D. = [-1270977, 427140, -1427543055, 579201840], r < 5 



a = (241530), N.D. = [502176, 175392, 1361800512, 14732928], r = 7 



a = (503412), N.D. = [11376, -3744, 745632, -292032], r = 5 



068 y2 + 300YX + 22500X2 - X^ + 894074112X - 19561912750080 

a = (014523), N.D. = [28656, -11808, -8478432, 2904768], r = 6 



069 y2 - 300YX + 22500X2 - X^ + 8940741 12X - 19561912750080 

o- = (230145), N.D. = [-30384, 11808, -6045408, 2904768], r = 6 



070 y2 - 60yX + 900X2 _ x'-'' + 48045312X - 189087436800 

a = (341250), N.D. = [-7344, 3744, -714528, 292032], r = 5 
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071 y^ - 960yX + 230400X2 - X^ + 13201887036672X - 17372103328571019264 
a = (502341), N.D. = [5909328, -1703520, 8039027520, -3454738560], r = 6 



072 Y^ + 600YX + OOOOOX^ - X^ + 4185188352X - 119899753943040 
a = (013452), N.D. = [52704, -12384, -22572864, 6390144], r < 7 



073 Y^ + 768YX + 147456^^ - X^ + 28465015272192X - 52273748736104702976 
a = (124503), N.D. = [4297392, -2075760, -8562479328, 3860913600], r < 6 



074 y2 - 4368yX + 4769856X2 - X^ + 102164509274112X - 521183246173070966784 
a = (235014), N.D. = [7567488, -3850560, 20709101952, -8825483520], r < 6 



075 y2 - 192yX + 9216X2 - X^ + 1067634432X - 20407326188544 
a = (340125), N.D. = [-32832, 12240, -6162912, 2496960], r = 4 

076 y2 + 120yX + 3600X2 „ js(^3 ^ 210573084672X - 33546735732363264 
a = (451230), N.D. = [-495072, 211680, -339888960, 124467840], r = 6 



077 y2 - 180yX + 8100X2 - X^ + 51432192X - 368371860480 
a = (125430), N.D. = [14448, -4752, -638496, 313632], r = 4 



a = (230541), N.D. = [169968, -66960, -19904832, 12454560], r = 5 

079 y2 + 396yX + 39204X2 - X^ + 388106595072X - 88686989876929536 
a = (341052), N.D. = [461712, -235440,459510624, -153977760], r = 5 



080 y2 - 60yX + 900x2 _ x'-^ + 48O45312X - 189087436800 

a = (214305), N.D. = [-7344, 3744, 745632, -292032], r = 5 



081 y2 - 3ooyx + 22500x2 - x^ + 894074112X - 19561912750080 

a = (325410), N.D. = [-30384, 11808, 8478432, -2904768], r = 6 



082 y2 + 300yX + 22500x2 - X-^ + 894074112X - 19561912750080 
a = (541032), N.D. = [28656, -11808, 6045408, -2904768], r = 6 

083 y2 - 60yX + 900X2 _x-^ + 48045312X - 189087436800 

a = (052143), N.D. = [11376, -3744, -714528, 292032], r = 5 

084 y2 - 3240yX + 2624400X2 - X^ + 2359518446592X - 1392039372309073920 
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a = (314025), N.D. = [502176, 175392, 1435465152, -14732928], r = 7 

085 Y^ - 1710yX + 731025X^ - X^ + 1041760268307X - 785868626594613294 
a = (425130), N.D. = [-1270977, 427140, 1468466145, -579201840], r < 5 

086 y^ - 396yX + 39204X2 - X^ + 16655166000X - 690091901769600 
<T = (530241), N.D. = [250140, -71280, 33275880, -18817920], r = 4 

087 Y^ - 612yX + 93636X2 - X^ + 52888734000X - 4252215735446400 
a = (041352), N.D. = [-243780, 110160, -126151560,44945280], r < 5 

088 Y^ + 4284YX + 4588164X2 - X^ + 252938243652912X - 2514799956670531410816 
a = (152403), N.D. = [23543292, -7847280,80180396856, -32581906560], r < 7 

089 y2 _ 5i84yx + 6718464X2 - X^ + 142783615579392X - 708085096419827994624 
a = (203514), N.D. = [13189008, -5387040, 28529516736, -13898563200], r = 4 

090 y2 - 60yX + 900X2 _x-^ + 767232X - 253808640 

a = (410325), N.D. = [-432, 288, 16416, -1728], r = 3 

:2 -x^- 

a = (032541), N.D. = [1008, -288, 7776, 1728], r = 3 



a = (143052), N.D. = [963, -360, -15444, 5400], r = 2 



a = (305214), N.D. = [-3348, 1440, 92448, -43200], r = 2 

094 y2 + 120YX + 3600X2 _ x'-^ + 81948672X - 327577374720 
a = (514320), N.D. = [-7776, 3168, -921024, 418176], r = 4 

095 y2 - 960yX + 230400X2 - X^ + 19321742592X - 1033886822584320 
a = (025431), N.D. = [91152, -6048, 38273472, 72576], r = 5 

096 y2 _ 4944yx + 6110784X2 - X^ + 14552162687232X - 21391150302252205056 
a = (130542), N.D. = [2867472, -340560, 5243236704, 44953920], r < 7 

097 y2 _ 3072yx + 2359296X2 - X^ + 2328753038487552X - 49786655684722942869504 
a = (241053), N.D. = [41436288, -19313280, -314701016832, 111012733440], r < 6 
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098 y^ - 2616yX + 1710864X^ - X^ + 3327107104512X - 3360005144504534016 
a = (352104), N.D. = [1503792, -403920, -2077894944, 906396480], r < 6 

099 Y'^ - 960yX + 230400X2 - X^ + 92277352334592X - 342160070800370356224 
a = (403215), N.D. = [-3525552, 3013920, -19974703680, 10813944960], r < 4 



100 Y'^ - 60YX + 900X2 - X^ + 767232X - 253808640 

a = (145230), N.D. = [-432, 288, 7776, 1728], r = 3 



As Bremner noticed in [S] , points in arithmetic progression seem to 
have a tendency to be independent. If we only take into account the 
ranks actually computed, we get an average over 4, when the average 
for random curves is known to be much smaller ([IDj). 

It should be noticed that example 039 is a remarkable case: no 
x-a.p. of length 6 were known in groups of rank one (see the last 
remark of 0), and |^ shows why this is so uncommon. The relation, 
if any, between the rank and the length of s.a.p. (or x-a.p.) seems a 
much harder problem to tackle. 
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